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I.  INTRODUCTION 


In  this  report,  we  consider  in  some  detail  various  approaches  to  the 
problem  of  processing  data  from  arrays.  This  problem,  which  is  frequently 
referred  to  as  the  combined  space -time  processing  problem,  will  be  discussed 
in  the  context  of  a  sonar  system.  The  general  results  are  applicable  to  array 
processing  in  any  area , 

In  Section  A, we  will  discuss  some  of  the  considerations  that  arise  in 
modelling  the  sonar  problem  and  indicate  some  of  the  implications  of  the  various 
models  .  In  Section  B,  we  develop  in  detail  the  quantitative  model  that  we  wul 
use.  In  Section  C,  we  outline  the  organization  of  the  report  and  summarize  the 
principal  results . 

Before  proceeding,  it  is  appropriate  to  mention  the  background  as 
sumed  of  the  reader.  A  knowledge  of  random  process  theory  at  the  level  of 
Davenport  and  Root^1)  is  necessary.  In  addition,  the  elements  of  statistical 
detection  theory  (e.g.,  Helstrom'2'  or  Ref.  1,  Chap.  14)  and  linear  filtering 
theory  (Ref.  1,  Chap.  11)  are  needed.  Matrix  notation  and  a  few  simple 
matrix  properties  are  used  (e.g.,  Hildebrand^)  0r  Bellman^4)).  In  various 
portions,  additional  background  it;  needed;  this  background  is  contained  in 
VanTrees<5). 


A.  PRELIMINARY  CONSIDERATIONS 


The  basic  system  of  interest  is  shown  in  Figure  1 .  The  waveforms 
are  received  by  an  array  of  hydrophones .  These  waveforms  contain  a  compo¬ 
nent  due  to  various  noise  sources,  which  may  be  within  the  hydrophones  or 
external  to  them.  In  the  case  of  active  sonars,  there  is  also  reverberation 
return.  If  a  target  is  present,  a  "signal”  component  will  be  added.  In  the 
active  sonar  case,  this  signal  is  a  reflection  of  the  transmitted  signal  from  the 
target.  In  the  passive  sonar  case,  the  signal  consists  of  sound  generated  by 
the  target  itself.  The  purpose  of  processing  it  to  obtain  information  about  the 
target  from  the  received  waveforms . 

The  first  problem  is  to  ievelop  a  suitable  model  for  the  signals  and 

noises . 
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FIGURE  1  SPACE-TIME  PROCESSING  MODEL 


The  received  waveforms  correspond  to  the  N  hydrophone  mrputs: 

r  (t).  r  ft).  -  •  -  r  ,(l)  It  ts  convenient  to  describe  these  as  a  single  vector 
1  2  N 

signal  r <t) 

I 

fj(t) 

r.,0) 


r(0  5 


rNW 


T  !  t  - T, 
1  1 


(Tj  is  the  initial  observation  time.  Tf  is  the  final  observation  time  ) 
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Similarly,  we  can  describe  the  signal  component  in  the  N  waveforms 
as  s(t)  and  the  noise  component  in  the  N  waveforms  as  n(t).  Thus,  when  a 
target  is  present, 


_r(t)  =  s(t)  +  n(t)  T  ^  t  iTf 

First,  we  consider  possible  methods  of  characterizing  the  signal 
component  s(t).  Four  possible  characterizations  come  to  mind: 

1 .  Deterministic  (Known)  Signals  -  In  an  active  system,  the  shape  of  the 
transmitted  signal  is  known.  the  channel  does  not  distort  the  signal 
and  the  target  acts  as  a  perfect  point  reflector,  then  the  shape  of  the 
signal  can  be  considered  known.  This  particular  model  is  rarely 
appropriate  for  the  sonar  case. 

2.  Deterministic  Signals  with  Unknown  Parameters  -  Normally  the  trans¬ 
mitted  signal  is  centered  around  center -frequency .  For  example,  the 
signal  illuminating  the  target  might  be 

s,  (t)  =  sin  jj  t 
I  c 

Even  a  simple  target  will  introduce  a  p^ase  angle  and  an  attenua¬ 
tion  which  will  be  unknown  to  the  receiver.  Thus,  the  reflected  signal 
might  be 


SR(0  .  VR  sin  (V  ♦  -R) 

where  V  and  fi  are  unknown.  Frequently,  it  is  reasonable  to 
R  R 

assume  that  ;  is  a  random  variable  with  a  unifuim  probability  density 
R 

over  the  interval  [0,2"  j.  The  constant  V  can  be  modelled  as  either 

R 

a  random  variable  or  as  an  unknown,  non -random  variable.  For  s 
simple  non -dispersive  channel  which  is  essentially  constant  ('  iring 
the  signalling  interval,  the  constants  can  include  its  attenuation  and 
phase  shift 


3.  Random  Signals  -  If  the  channel  or  target  fluctuates  while  being  illumi¬ 
nated  by  the  signal,  then  the  signal  will  be  distorted.  Since  these 
fluctuations  arc  inherently  •"andom,  a  convenient  approach  is  to  view 
the  signal  as  a  sample  function  from  a  random  process. 
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In  the  passive  sonur  case,  the  signal  emitted  from  the  tar¬ 
get  is  caused  by  a  variety  of  things  such  as  engine  noise  or  propeller 
noise.  Here,  the  signal  source  is  inherently  random. 

In  both  of  these  cases,  one  usually  knows  enough  about  the 
physical  situation  to  be  able  to  specify  the  approximate  second -moment 
characteristics  of  the  process  (i.e.,  the  mean  value  function  and  the 
covariance  function).  Frequently  the  physical  origins  of  the  random¬ 
ness  are  such  that  one  can  assume  the  process  is  a  Gaussian  random 
process.  In  this  case,  the  second -moment  properties  provide  a  com¬ 
plete  characterization. 

Once  we  have  assigned  statistical  properties  to  the  signal, 
we  would  expect  that  any  optimum  processing  schemes  we  design  will 
be  good  on  the  average.  In  other  words,  any  particular  time  the 
experiment  is  conducted  the  performance  may  be  good  or  bad.  but 
when  averaged  over  the  assumed  signal  ensemble  the  performance  will 
be  optimum . 

In  many  cases  it  is  difficult  to  assign  statistical  properties 
to  the  signal .  Alternately,  we  may  want  to  design  a  test  that  will  be 
optimum  for  the  particular  signal  present  (as  contrasted  to  a  signal 
ensemble).  In  nese  situations,  a  fourth  model  is  appropriate. 

4.  Non -Random,  Unknown  Signals  -  In  this  case,  we  assume  that  the  N 
signal  inputs  s.(t),  s  2  (t) ,  •••  s^  (t)  are  unknown.  We  want  to  design 

processors  that  are  good  for  the  particular  set  that  is  present  (not 
some  average  signal  set). 

Next,  we  consider  characterizations  for  the  noise.  As  mentioned 
earlier,  the  possible  sources  include  I)  hydrophone  noise,  2)  ambient  noise, 
and  3)  reverberation  noise  in  active  systems.  The  first  two  arise  lrom  the 
combined  effect  of  many  small,  independent  sources  .  According  to  the  Central 
Limit  Theorem,  one  can  model  these  as  Gaussian  (normally  stationary)  random 
processes . 

Reverberation  noise  is  caused  by  the  reflection  of  the  transmitted 
signal  from  various  objects  in  the  ocean.  Its  properties  will  be  a  function  of 
the  transmitted  signal  shape  and  the  reflection  mechanism.  A  detailed  deri¬ 
vation  of  a  possible  mathematical  model  is  given  by  Kelly  and  Leiner^)  (see 
also  Refs.  7  and  8).  If  one  assumes  a  large  number  of  small  reflectors,  one 
is  led  to  a  Gaussian  process  model  which  is  non -stationary  and  has  character¬ 
istics  that  depend  on  the  transmitted  signal. 
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Accordingly,  we  will  assume  that  th°  noise  is  a  sample  function  from 
a  Gaussian  random  process  with  known  statistical  properties. 

In  addition  to  developing  a  suitable  model  for  the  signal  and  noise 
environment,  we  must  establish  the  goal  of  our  processing  system  and  specify 
a  criterion  to  measure  how  closely  the  system  achieves  this  goal. 

We  will  discuss  some  possible  criteria  for  the  various  signal  models 
suggested  above. 

1 .  Case  1 

When  the  shape  of  the  transmitted  signal  is  known,  the  question  of 
interest  is  normally  whether  or  not  the  target  is  present.  This  is  a  binary 
hypothesis  testing  problem  .  The  received  waveform  under  the  two  hypotheses  is: 

H,  :  r  (t)  -  s  (t)  +  n  (t)  T.  £t  T,  (target  present) 

Hr  :  r(t)  -  n  (t)  T.  M  £  T#  (target  absent) 

u  —  —  it 


Using  either  a  Bayes  or  Neyman -Pearson  criterion,one  is  led  to  a 
likelihood  ratio  test.  One  operates  on  the  waveform  r(t)  to  construct  the  function, 


A(r(t))  ^ 


"P  (£  (t) |  Hj)" 

"p  <£(t)|H0r 


(The  quotation  marks  emphasize  that  we  must  be  careful  about  the 
meaning  of  these  expressions.) 

One  then  compares  A  to  some  threshold  r,  which  is  chosen  as  a 
function  of  cost  in  the  Bayes  case  or  to  achieve  a  desired  Pp  in  the  Neyman  - 
Pearson  case . 

The  resulting  processor  for  known  s_(t)  turns  out  to  be  a  matrix 
linear  filter  with  N  inputs  and  1  output. 

Alternately,  one  cannot  assume  the  noise  is  Gaussian  if  only  the 
second -moment  properties  of  the  noise  are  known.  Then,  one  does  not  have 
enough  information  to  construct  A(r(t)).  A  plausible  approach  in  this  ase  is 
to  require  a  linear  processor  (once  again,  with  N  inputs  and  1  output)  and  try 
to  maximize  the  output  signal -to -noise  ratio, 
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2 

(output  due  to  s  (t)  at  T,) 

S  L  _ - _ 

N  E  £ (output  due  to  n  (t)  at  Tf)  J 

(The  symbol  E  denotes  expectation.) 

It  is  straightforward  to  show  that  these  two  criteria  lead  to  identical 

results . 


2 .  Case  2 


If  the  unknown  parameters  are  random  and  the  problem  is  one  of 
detection,  the  procedure  is  identical  to  Case  ! .  However,  if  the  parameters 
are  non -random  or  if  we  want  to  estimate  them,  the  procedure  must  be  modified. 


.3 .  Case  3 

For  random  signals,  two  possibilities  exist. 

First,  if  detection  is  the  problem  of  interest,  we  are  once  again  led 
to  a  likelihood  ratio  test.  To  construct  the  likelihood  function,  it  is  necessary 
to  characterize  the  signal  process  completely.  The  most  common  process 
model  is  a  Gaussian  model.  As  discussed  above,  the  second -moment  charac¬ 
terization  then  provides  a  complete  description.  In  our  detailed  discussion,  we 
will  encounter  certain  cases  in  which  the  Gaussian  signal  assumption  will  be 
valid.  The  resulting  processor  is  a  quadratic  device. 

Second,  it  is  often  desirable  to  estimate  what  the  signal  component 
of  the  input  is.  Since  both  the  signal  and  noise  are  sample  functions  of  a  ran¬ 
dom  process,  a  minimum  mean -square  error  estimate  is  appropriate.  We 
denote  this  MMSE  estimate  by  s  (t).  Under  the  assumption  that  the  signal  and 
noise  process  are  both  Gaussian,  s(t)  is  obtained  by  using  a  linear  processor. 
On  the  other  hand,  the  Gaussian  signal  process  assumption  is  not  invoked,  the 
form  of  the  processor  must  be  specified.  If  we  ask  for  the  best  linear  MMSE 
estimate,  we  can  solve  the  problem  using  only  second -moment  properties  of 
the  process .  The  resulting  estimates  are  the  same  in  both  cases . 
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4.  Case  4 


For  non -random  but  unknown  signals,  the  desired  procedure  is  less 
obvious.  There  are  two  possibilities,  which  are  analogous  to  those  outlined  in 
Case  3. 


One  can  construct  a  generalized  likelihood  ratio  test: 


max 

A(£(t))  L  s  (t) 
6 


*’P  Lr  (0  |  Hj  ]" 


”P  [r(t)  |  HQ  ]" 


and  compare  A„  to  some  threshold  adjusted  to  give  the  desired  false  alarm 
probability.  Tne  numerator  is  found  by  making  a  maximum  likelihood  estimate 
of  s  (t)  and  substituting  it  into  the  probability  density.  This  test  has  no  claim 
to  optimality  but  is  intuitively  logical  and  frequently  performs  well . 


The  analogy  to  a  minimum  mean -square  error  estimate  of  a  sample 
function  is  the  maximum  likelihood  estimate  of  a  non -random  function.  If  one 
does  not  invoke  the  Gaussian  assumption  on  the  noise,  then  it  is  appropriate  to 
ask  what  sort  of  estimate  of  s  (t)  can  be  obtained  using  a  linear  filter.  An 
appropriate  criterion  might  be  to  require  the  output  to  equal  s  (t)  exactly  if  the 
noise  were  absent.  Subject  to  this  constraint,  we  design  the  linear  filter  to 
minimize  the  distortion  due  to  noise.  The  resulting  estimate  is  found  to  be  the 
same  as  the  maximum -likelihood  estimate. 


In  this  section,  we  have  discussed  some  of  the  considerations  that 
are  involved  in  choosing  an  appropriate  mathematical  model  for  the  physical 
problem  of  interest.  We  will  concentrate  our  attention  on  Cases  3  and  4;  the 
noise  will  be  modelled  as  a  sample  function  from  a  Gaussian  random  process, 
and  the  signal  will  be  modelled  either  as  a  sample  function  from  a  Gaussian 
random  process  or  as  a  non -random,  but  unknown,  waveform. 

We  shall  now  specify  this  model  in  detail. 


B.  QUANTITATIVE  MODEL 

The  general  model  is  easily  stated.  The  received  waveforms  of 
interest  are  r ^  ( t) ,  •  •  ■  r^  (t).  First,  consider  the  noise  components, 
nj  (t),  ' ' '  n n  ( t) .  We  denote  these  by  the  vector  n(t).  The  following  proper • 
ties  are  assumed: 

1 .  Each  noise  nj  (t)  ;  i  =  1,2,  •  •  •  N  is  a  sample  function 
from  a  zero -mean  Gaussian  process . 
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2 .  The  noise  processes  are  jointly  Gaussian . 

3.  Statements  (1)  and  (2)  are  equivalent  to  saying 
that  n  (t)  is  a  zero -mean  vector  Gaussian  process  . 


(Observe  that  the  covariance  function  and  correlation  function  are  the  same 
because  of  the  zero -mean  assumption.) 


5.  We  further  assume  that  each  noise  function  contains 
a  non -zero  "white"  component  which  is  independent 
of  the  remaining  noise  function  and  of  the  "white” 
components  in  the  other  noise  waveforms  .  Thus , 
we  may  write 

n  (t)  =  w  (t)  +  n  (t) 

i  i  c} 

where  w^(t)  represents  the  white  component  and 
n^  (t)  is  the  remaining  colored  noise  component. 


K  (t ,  u)  =  K  (t ,  u)  +  K  (t ,  u) 
n,  w,'  n 


6  (t-u)  +  Kn  (t . u) 
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For  algebraic  simplicity,  we  assume  that  the  white 
components  in  each  process  are  equal.  (The  general 
case  is  a  trivial  modification.) 

The  covariance  matrix  of  the  vector  process  is: 

K  n(t>  u)  =  y  fc(t-u)l+  Kc(t,u) 

where  J.  is  the  identity  matrix  and  K  (t,u)  is  the 
covariance  matrix  of  the  colored  noises  . 


Physically,  this  white  component  may  correspond  to 
internal  noise  in  the  hydrophone  and  its  associated 
circuitry.  Since  its  bandwidth  is  wider  than  any 
signals  of  interest,  we  may  model  it  as  a  process  with 
a  flat  spectrum  at  all  frequencies  . 


We  assume  that  the  colored  noise  component  has  finite 
mplies  that 

K  (t,u)  KT(t,u)dtdu  <  « 
c  c 


Observe  that  we  do  not  need  to  be  explicit  about  the 
noise  sources  .  Any  Gaussian  noise  source  which  leads 
to  the  same  vector  covariance  function  will  be  treated 
alike.  Later,  in  some  examples,  we  will  see  how 
various  noise  fields  give  rise  to  particular  covariance 
functions . 

The  signal  is  also  characterized  by  a  vector  s  (t). 

Under  the  Gaussian  assumption,  we  denote  'ts  covariance  by 
Ks  (t,u)  and  assume  it  is  zero -mean .  It  does  not  contain  a 
white  component.  For  simplicity,  we  assume  it  is  indepen¬ 
dent  of  the  noise  process  . 

Observe  that  this  would  not  be  true  if  the  channel  or  target 
fluctuated  while  being  illuminated  by  the  signal,  because 
the  transmitted  signal  would  influence  both  the  returned 
signal  process  and  the  reverberation  noise.  The  modifi¬ 
cation  to  include  this  coupling  is  straightforward. 


energy.  This 
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For  the  general  developments,  no  further  restrictions  on  the 
components  of  s(t)  are  needed.  A  special  case  that  we  will  emphasize  in 
most  of  our  examples  is  one  in  which  each  signal  sj(t)  is  a  shifted  version 
of  the  same  signal . 


s.(t)  =  s(t  -  t.)  i  =  1,2,  •  •  •  N  (1) 

where  the  t  .  are  known. 

1 

Physically,  the  simplest  case  this  might  correspond  to  is  a  point 
target  and  plane -wave  propagation.  Neither  of  these  assumptions  are 
necessary. 

Under  the  non -random  signal  assumption,  s(t)  is  simply  an 
unknown  vector.  Once  again,  we  will  emphasize  the  case  described  by 
Eq.(l). 


C .  ORGANIZATION  AND  PRINCIPAL  RESULTS 


In  Chapters  II  and  III  we  study  the  random  signal  case  in  some  detail. 
First,  we  find  the  form  of  the  optimum  processor  for  the  general  case.  Then 
we  look  at  some  interesting  special  cases  and  evaluate  their  performance. 

In  Chapter  IV  we  study  the  non-random,  but  unknown,  signal  problem. 

In  Chapter  V  we  look  at  the  special  case  in  which  the  signal  vector 
satisfied  Eq.  (l)and  the  noise  is  homogeneous  (i.e. ,  Kn.(t,u)  is  not  a  function 
of  i).  The  notion  of  array  gain  is  encountered  and  its  s\gnificance  discussed. 

In  Chapters  VI  and  VII  we  consider  some  particular  examples  of  dis¬ 
tributed  and  directional  noise  fields  . 

Finally,  in  Chapter  VIII  we  summarize  the  results  and  suggest  some 
future  work . 

The  principal  result  is  the  demonstration  that,  for  a  large  class  of 
problems,  the  basic  structure  of  the  optimum  processor  does  not  depend  on  the 
signal  case  of  interest  nor  on  the  criterion  used.  Specifically,  we  will  arrive 
at  a  receiver  in  which  the  only  matrix  operation  is  invariant  to  the  above  assump¬ 
tions  and  the  solutions  to  all  of  the  above  problems  appear  at  various  points  of 
the  receiver. 

A  secondary  result  is  some  insight  into  situations  in  which  optimum 
processing  may  be  worthwhile. 
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II.  GAUSSIAN  SIGNALS  IN  GAUSSIAN  NOISE 


In  this  chapter  we  derive  the  optimum  detector  for  the  random  Signal 
process  described  in  Case  3. 

The  basic  derivation  is  not  restricted  to  the  array  processing  case. 
Historically,  the  essential  results  were  first  obtained  by  Price  in  1954  during 
his  studies  on  multiple  scatter  links  for  communication ^  In  1959  the 
problem  was  studied  in  a  different  context  bvWolf.^)  Recently,  Bryn^)  re¬ 
derived  the  results  for  the  special  case  of  stationary  processes  and  infinite 
time  intervals.  Middleton  and  Groginsky(^  have  also  studied  the  array 
problem . 


A.  DERIVATION  OF  LRT 


As  discussed  in  the  introduction,  the  target  or  channel  changes 
appreciably  during  the  observation  interval  in  many  active  sonar  situations  . 
Thus,  even  though  the  transmitted  signal  is  completely  known,  the  effect  of  the 
transmission  path  or  the  target  reflection  mechanism  causes  the  returned  wave 
form  to  have  a  random  behavior.  With  passive  sonar,  the  actual  generated 
waveform  has  a  non -deterministic  character. 

A  Gaussian  model  is  suitable  for  this  situation.  In  the  absence  of  a 
target,  the  received  waveform  will  consist  of  contributions  due  to  the  various 
types  of  noises  discussed  above.  As  before,  we  call  this  hypothesis  H^: 

£(t)  -  n(t)  [Hq]  (2) 


We  will  assume  that  there  is  a  non-zero  amount  of  noise  generated  at 
each  hydrophone.  Tins  noise  is  statistically  independent  at  the  various  elements 
in  the  array  and  is  assumed  to  be  "white”  over  the  frequency  ranges  of  interest. 
To  emphasize  this  white  noise  component,  we  will  occasionally  write 


n(t)  =  w(t)  +  n  c(t) 


(3) 
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The  covariance  matrix  under  is: 


KQ(t,u)  =  E  |^r(t)rT(u):  HQ  is  true j 


=  E  [n(t)  nT(u)  ]  =  Kn(t,u) 


(4) 


AN0 

=  I  +  K  (t.u) 
z  0 


Under  hypothesis  Hj,  there  is  an  additional  random  component  due 
to  the  target: 


1  (t)  =  s  (t:  f.)  +  n(t) 


(5) 


The  signal  s  (t:  f<)  is  a  sample  function  from  a  Gaussian  random 
process  .  For  our  present  purposes,  it  is  adequate  to  assume  that  the  process 
is  zero -mean. 

The  covariance  under  Hj  is: 


K  j  (t.u)  =  E  ^  [  r(t)  -  m  j(t) ]  i.rT (u)  m  ^(u)  ]  j 


(6) 


£  K  (t,u)  +  K  (t , u) , 
—  s  —  n 


since  we  assumed  the  mean  is  zero. 

To  solve  the  det<“Ction  problem,  we  must  compute  the  likelihood  ratio 
and  compare  it  to  a  threshold.  The  likelihood  ratio  test  is 


•’>  ( r  <t))  = 


P  [ £ (t)  Hj  ] 
p  L  £  (t)  H0  ] 


(7a) 
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There  are  several  ways  of  formulating  this  likelihood  ratio.  A 
method  chat  leads  to  simple  interpretation  of  the  result  is  to  construct  A  (r  (t)) 
as  a  ratio  of  two  fractions: 


A(r(t)) 


P  [l(t)  I  Hl  ] 

p  [r(t)  |  w(t)  only  J 


p  b>  I  H0  ] 

P  [r(t)  I  w(t)  only  j 


(7b) 


It  is  easy  to  show  that  both  terms  exist.  Since  they  have  a  similar 
structure,  it  is  adequate  to  study  the  first  term  in  detail. 


[r(t)  1  Hj  ] 


p  jjr(t)  ,  w(t)  only  J 


lim 

N 


1  N  /  T( 

(2*)n/2  TT  «p(-  ?  11 * 

1  =  1  \  T, 


du  r  (t)Qj  (r.u)  r(u) 


N 


(2n) 


ntt  TT  x>xpIn7  l/(t) 


i  =  1 


I  r(t)  dt^ 


(8) 


where  the  are  »he  scalar  eigenvalues  in  the  vector  Karhunen*Lo£ve 
expansion^  and  Q  y  (t,u)  is  the  inverse  matrix  kernel. 


The  inverse  matrix  kernel  satisfies  the  equation 

Tf 

J  K  j(t,u)  Qj(u,z)  du  •  *(t-z)J_  <9) 

Ti  " 


In  terms  of  the  vector  eigenfunctions 

Nr 

.  ( 

i  =  l 


*  ,  N0  c  ,  n  T 

— , <f,u)  =  i  (t  +  \  V  (u) 


(10) 


or 


2  +  i 


(ID 


t  The  vector  Karhunen -Lo£vc  expansion  was  developed  by  Kelly  and  Root. 


(15) 
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Taking  the  .ogarithm  of  Eq.  (8).  we  have: 


X  *  1  . 

In  A(r(t))  =  +  j  Cl  ^  ff  vV0  *  \  jj  dt  du  IT(t)  [O.1  (t,u)  '  fT-  f  jl1 
i=  1  t)  7  o 

1  (12a) 

The  first  term  represents  a  bias  which  is  not  dependent  on  the  received  wave* 
form;  the  second  term  represents  the  operation  on  the  received  waveform. 

We  will  first  develop  two  canc  nic  forms  for  the  operation  on  the  re¬ 
ceived  data  and  then  develop  a  convenient  expression  for  the  bias  term .  We  will 
denote  the  period  term  by  the  symbol  L.  The  total  test  will  consist  of  ;wo  terms 
like  Eq.  (12a). 

Using  Eqs.  (7b)  and  (12a),  we  see  that  the  likelihood  ratio  test  becomes: 

T( 

-  j  J  dt  du  rT(t)  j (t . u)  -  |  ;  (t-u)J  r  (u) 

Tj  ° 

Tf 

+  j  J  dt  du  rT  (t)  [(^(t.u)  -  ~  ^  Mt-u)^  r(u) 


(12b) 


1  *  ,  2  c,\  1  *  r  2 

<  r’  2  L  .,+  NTi>2  N"  '?J 


We  will  initially  investigate  the  first  term  on  the  left-hand  side  of 


Eq.  (12b). 


1 .  Receiver  Structure  «1 


To  find  a  convenient  form  wc  divide  the  covariance  function  into  two 


parts: 


K,  (t  u)  =  —  •  (t  u)  I  *  K  (t.u) 

-I  2  -  -c, 
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We  want  to  find  a  solution  to  Eq.(9)  with  the  following  form: 


Q,  (t.u)  =  ~  {  Mt-u)l  -  hc  (t,u)  ]• 


Substituting  into  Eq .  (9),  we  have 


0  f 

K  (t.u)  =  +  h  (t.u)  +  K  (t.z)  h  (z,u)  dz  T  s  t,u  s  T 
— c ,  z  —c ,  J  —  c ,  —  e ,  i 


A  formal  solution  can  be  obtained  by  writing 


or 

h  (t.u)  =7  h;  V  (t)  V  T  (u) 
-c.  Z-.  i  i  l 


Substituting  into  Eq.(15)  and  using  Eq .  (10),  we  have: 


2  .  c  . 

hl  _  No  ‘‘ 

'  l  +  ir'Cl 

nj  1 


Using  Eqs .  (12)  and  (14).  we  see  that  the  first  quadratic  form  on  the 
left  side  of  Eq .  (12b)  is: 

Tf 

+  J_  jj  df  dtr  rT(t)  h  (t.u)  r  (u)  (|.S) 

No  ’  "CI 

« i 

Now  the  inner  integral  in  Eq .  (18)  >s  famihar  from  optimum  filter 
theory  (eg..  Ref.  Chap,  h) . 


|  h  (t.u)  r  (u)  do  =  s  (t  :  H  •)  a  n  (t  :  H - ) 

1  —  c .  —  —  —  c 
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Thus,  the  operation  represented  by  the  inner  integral  in  Eq .  (18)  is 
equivalent  to  an  estimation  of  the  non -white  portion  of  the  input,  assuming  that 
Hj  is  true. 

Proceeding  in  an  identical  manner  with  the  second  term  in  Eq .  (12b), 
we  obtain  the  receiver  shown  in  Figure  2. 


FIGURE  2  RECEIVER  STRUCTURE  *1 


We  observe  that  there  is  only  one  filter  in  the  receiver. 

h  .  (Uu)  _A  [h  (t.u)  -  h  (t,u)  1  (20) 

-A  -  Nq  L-Cj  c o  J 

One  can  also  show  that  this  filter  satisfies  the  following  integral  equation 
(see  App.  A  for  proof): 

Tf 

+  K  (t.x)  =  f  K.  (t.u)  h  ,(u,z)  Kn  (z,x)  dudz  T.  st,  x  £T,  (21) 

Ti 

or,  using  Eq  .  (6), 

Ks  (t.x)  =  jj  [jCs(t.u)  +  K  n  (t .  u)  J  h  ^(u ,  z)  [k  n  (z ,  x)  j  du  dz  (22) 
Ti 
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2  .  Receiver  Structure  #2 


obtained. 


A  second  interesting  interpretation  of  the  optimum  receiver  can  be 
We  divide  h  (u,z)  into  two  parts: 


h,(u.z)=  g(u,w)  K  *  (w ,  z)  dw 
—  a  «l  —  —  n 


(23a) 


T; 


where 


l 

s 

T. 


(t.u)  (u,z)  du  =  6(t-z)| 


(23b) 


Then,  substituting  into  Eq.  (22),  we  obtain  an  integral  equation  that 
g(u,z)  must  satisfy: 


-s  (t,X)  =  I  [-s(t,U)  +  -n  (t,u)  ]l(u’x)  du 


T.  - 1 ,  x  £  Tf 


(24) 


This  equation  is  familiar;  the  matrix  g(u,x)  is  related  to  the  matrix 
filter  one  uses  to  find  the  minimum  mean -square  estimate  of  s(t),  given  the  input 
r(t).  Specifically, 


Tf 

£(u)  =  J  gT(x,u)  r  (x)  dx  (25) 

Ti 

The  test  statistic  is: 

Tf 

L  =  JJ  dt  du  rT (t)  h  ^(t.u)  r(u)  (26) 
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(The  factor  of  1/2  is  absorbed  in  the  threshold  .) 


Substituting  Eq.(23)  into  Eq.  (26), we  have: 

Tf 

L  =  I  dt  du  dv  (t)  g  (t,v)  K  1  (v,u)  r  (u)  (271 

c’  J  w  —  —  —  — . 

Ti 

W  riting 

Tt  . 

z  (v)  =  f  (v,u)r(u)du  (28) 

Ti 

we  have 

Tf 

L  *  j  £  ^  (v)  £  (v)  dv  (29) 

Ti 


The  receiver  shown  in  Figure  3^  is  an  alternate  version  of  the 
estimator -correlator  of  Figure  2  that  we  find  useful  in  the  sequel . 


FIGURE  3  ESTIMATOR  CORRELATOR  RECEIVER: 

FORM  #2  OF  THE  OPTIMUM  RECEIVER 


T 

t  The  symbol  ©  denotes  a  dot  product  A  B  .  The  double  lines  indicate  vector  signals  . 
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The  vector  z(t)  also  has  a  familiar  interpretation  from  the  active 
sonar  case  (e.g.,  Van  Trees, Ref.  16  or  5,  Chap.  4).  It  is  the  input  to  the 
correlator,  as  shown  in  Figure  4.  The  other  input  is  f(t),  the  known  signal. 
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Then 


(z,  u)  r  (u)  du 


(31) 


If  we  define 


Tf 

x  (z)  a  J  k"^  (z,u)  r  (u)  du  , 
~  Ti 


(32) 


then  the  two  inner  integrals  correspond  to  the  squared  magnitude  of  x(z). 


Tj 


f 

Ti 


L  =  J  dz  [  |  J  rT  (t)  k  (t, z)  dt  |2J 


(33) 


The  resulting  receiver  is  shown  in  Figure  5. 


FIGURE  5  FILTER-SQUARER  RECEIVER: 

FORM  #3  OF  THE  OPTIMUM  RECEIVER 
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The  two  remaining  problems  with  the  receiver  are: 

1 .  A  convenient  closed -form  expression  lor  the 
infinite  series  on  the  right  hand  side  of  Eq .  (12b) 
that  is  part  of  the  bias.  (This  term  is  evaluated 
in  Appendix  B.) 

2.  Some  measure  of  performance  o'  the  receiver. 
This  problem  is  discussed  in  the  next  section. 


B.  PERFORMANCE  OF  OPTIMUM  RECEIVER 


The  performance  is  difficult  tc  compute  in  the  general  case.  A 
quantity  that  provides  a  reasonably  good  indication  of  performance,  particularly 
when  the  input  signal -to-noise  level  is  low,  is: 


H 

f  "  1 

X 

J 

-  E 

-L  1 

H0]}2 

Var 

[LI 

I  1 

o 

X 

(34) 


This  corresponds  to  an  output  signal -to -noise  ratio.  Its  value  as  a  performance 
indicator  has  been  discussed  in  detail  by  Priced17) 

The  calculation  is  straightforward: 


L  = 


r  (t)  h  (t,u)  r  (u)  dt  du 


(35) 


4 


E  l|  h1]=e1  {Jj  £.T(t)  il^  <l'u)  £<u)  dt  du  } 


(36) 


But 


So 


AT  B  A  =  Tr  [  B  (A  A  T)  ] 


E  [l  |  Hj  ]  =  Tr  jj h^(t.u)  Kj  (t,u)  dtdu 


(37) 


(38) 
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Similarly, 


l 

E  [l  i  H()  j  =  Tr  Jj  h  ^(t,u)  K ^  (t.u)  dt  du 


(39) 


Thus,  the  term  in  the  numerator  of  Eq.  (34)  is: 


'  „  4  E  [L  '  H,  ]  -  E  [l  I  H0  j 


(40) 


Tf 


=  Tr  jj  Ks  (MO  dt  du 


(41) 


The  mean -square  value  follows  in  the  same  manner: 


Tf  Tf 


E  L2  ,  Hq  =  Eq  j  ■  •  J  rT(t)  h  ^(t.u)  r  (u)  rT  (x)  h  jx.y)  r  (y)  dtdudxdy 


(42) 


Using  the  factoring  properties  of  Gaust..~n  variables,  it  follows  that 


ell 


H 


0  J 


■  H1 


+  2  E 


0 


Tf  Tf 

j  •  •  j  rj  (t)  h  .(t.u) 


(u,x)h  .(x,v)  r(v)dtdudxdy 

'  (43) 


=  tE  lL  H„  ]  I2 


Tf 

I  h  (t.u)  K  (u.x)  h  .(\.v)  K  (v.t)dt dudxdv 

w  —  L  —  U  —  L.  —  U 

T, 


The  variance  is  just  the  last  term. 
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Therefore 


d 


2 


TfTf  *|2 

Tr  j  j  h  (t.u)Ks  (t.u)dt  du 
T  T 

_  1  i _ _ _ 

Tf  Tf 

2  Tr  J  •  *  J  h  L  (t,u)  K0  (u,x)  h  (x,y)  Ky  (y.t)dtdudxdy 
T,  T, 


(44) 


The  denominator  can  be  simplified  by  using  the  expression  in  Eq.  (23) 
for  h  ,  (t,u)  and  h  (x,y)  and  observing  that 

— 1~  —  a 


(v.u) 


Kq  (u,x)  du  =  \ 


4<V-X) 


(45) 


Performing  the  integrals,  we  obtain: 


d  =  2  Tr 


i 


g(t,x)  g  (x , t)  dx  dt 


T, 


(46) 


C .  COMMENTS 


In  this  chapter  we  derived  the  structure  of  the  optimum  receiver  for 
the  case  of  Gaussian  signals  in  Gaussian  noise.  To  implement  it.  we  must 
solve  Eq.  (22)  or  (24),  which  follow. 


Tf  _ 

(t,\)  =  jj  L-s  <l'U>  +  -n  (t,u)  I  (u.z)  |_5n^z,x^  du  T.  st.x'-Tf  (22) 

T, 
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or 


Tf 

K  (t.x)  =  J  i  K  .  (t.u)  +  K  (t.u)  ]  g  (u , x)  du  T.  *  t,  x  •  T,  (24) 

»  v  i  ”  s  *—  n  j  —  t  f 

Ti 

We  also  found  that  the  output  signal  -to  -noise  ratio 

Tf 

d2  =  2  Tr  jj  g(t.x)  g(x,t)  dx  dt  (46) 

Ti 

depended  on  the  solution  to  Eq.  (24).  For  arbitrary  observation  intervals  and 
random  process  statistics  the  solution  is  difficult.  Explicit  solutions  can  be 
obtained  in  several  cases  of  importance: 

1  .  If  the  signal  and  noise  processes  are  stationary 
with  rational  spectra,  one  can  transform  the 
integral  equation  into  a  differential  equation, 
solve  it,  and  substitute  the  solution  into  the 
integral  equation  to  satisfy  the  boundary  condi¬ 
tions.  This  case  is  conceptually  straightforward 
but  extremely  tedious. 


2.  When  the  noise  is  "white,’'  K  (t,  u)  can  be  written  as: 

—  n 

N0 

K  (t.u)  =  —  I 
—  n  l  — 

If  the  largest  eigenvalue,  >  ,  of  K  (t.u)  is  much  less 
than  N()/2,  a  solution  follows  easilv.  This  is  commonly 
referred  to  as  the  '  threshold''  or  "coherently  undetec¬ 
table"  case.  One  can  show  that 

h  (t.u)  >  K  (t.u). 

1.  If.  m  addition  to  the  conditions  of  case  (1),  we  assume  that 
T  -  T  -  T 

f  i 

is  long,  then  a  verv  simple  solution  can  be  obtained 
using  Fourier  transform  techniques.  Tins  case  is 
appropriate  in  most  siwiar  problems  and  is  developed 
in  detail  in  the  next  chapter. 
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III.  LONG  OBSERVATION  INTERVALS.  STATIONARY  PROCESSES 


A.  GENERAL  CASE 


In  this  case,  simple  solutions  for  the  integral  equations  can  be  ob¬ 
tained.  We  assume  that  the  signal  and  noise  processes  are  wide-sense  station¬ 
ary  and  that  the  observation  interval  is  long. 

First,  we  solve  the  equation  for  the  optimum  processor  filter 
hA(u,z).  We  let  Tj  =  -  T/2  and  Tr  =  +T/2.  Then,  by  Eq.  (21). 

+T/2 

K  (t.x)  =  >1  K,(t,u)hA(u.z)  K  (z.x)du  dz  (47) 

— s  j  J  —  l  —  c  —n 

-T/2 

We  then  multiply  by  (x.  y)  and  integrate  with  respect  to  x.  In 
addition,  we  can  write  the  covariance  functions  in  terms  of  differences  of  their 
arguments: 


+T/2 


+T/2 


-l 


(t  -  x)  K 

>  -n 


(x  -  y)  dx  = 


Kj(t  "  u)  h^(u  -  y)  du  (48) 


-T/2 


-T/2 


Now  we  define  the  Fourier  transforms  of  the  various  matrices.  For 

example. 


S  (x)  -  '  dVj*!?IC  (t) 
— n  «  -n 


(49a) 


and 


dx 

2” 


S  (x) 


(49b) 
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Observe  that 


sV)  =  sV> 

-n  - n 


(4“c) 


- 1 

Next  we  write  K  (')  and  h.(-)  in  terms  of  their  transforms. 
—  n  — i- 


+  T/2 

dx  K  (t 

—  S 

-T/2 


+  ® 


-n 


(SO) 


+  T/2 


K  ^(t  -  x)  dx 


+  jx(x  -  y)  tlx 


Or- 


-T/2 


Re^roupin^  terms. 


di  +jx<t  v) 


>- 


+  T/2 


-  T/2 


K  (t  -  x)  e 

-  s 


Ji<t  -  x) 


dx 


S 

-n 


(0 


(SI) 


dx 

2- 


o(t-y) 


+  T  2 


T  2 


K^t-xjc 


-J  (t  x) 


dx  HU) 


As  T  -  the  terms  in  the  parentheses  approach  Sx(i)  and  SjU) 
respect  t\elv .  Therefore. 


S  (i)S  "(»>  -  S  ( t)  H  .U) 

— s  — n  - 1  —  w 


02) 
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or 


H  ,(x)  =  S,'l('ii)S  (x)  S  (x) 
—  -I  — s  -n 

From  Eq.  (23).  it  follows  directly  that 


(S3) 


G(x-)  =  S  j  1  (jj)  S^(ji) 


(54a) 


and 


GT(x)  =  S*(x)  S  1  *(x) 

_  — s  —l 

This  result  is  well-known  from  unrealizable  Wiener  filter  theory.  The  output  is: 

l s(x)  =  GT  Ox)  R(jx)  (54b/ 


The  test  statistic  is: 

+T/2  +T,  2 

L  =  I  r  (t)  h  ,(t  -  u)  r(u)  dt  du 

W  J  —  — u  — 

-  T/2  -T/2 


This  can  be  expressed  in  the  frequency  domain  as: 


(55) 


T 

L  =  '  R*7(x)  H.(x)  R(x)  ~ 

*  — c  —  L 


(56) 
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tHcrc  R(x)  is  the  Fourier  transform  of  r(t).  To  be  correct,  one  should  use  the 
integrated  transform  but  the  ordinary  transform  is  more  familiar  so  we  use  it. 
A  script  £s(x)  denies  the  F  >uricr  transform  of  s(t).  while  an  ordinary  S^fx) 
denotes  the  power  density  spectrum . 
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Two  receiver  configurations  are  shown  in  Figure  6. 
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B.  SPECIAL  CASES 


1 .  Threshold  Case 

In  the  threshold  case  the  input  signal -to-noise  ratio  is  low.  If 

S  ..(id)  <<  S  ..(uu)  (58) 

s,  ij  n,  lj 

for  all  i,  j  and  uu,  then  we  call  the  problem  a  threshold  problem .  (Eq.  (58)  is  a 
stronger  condition  than  necessary;  some  weaker  conditions  will  arise  later.) 
This  equation  simply  says  that  each  element  in  the  signal  matrix  is  smaller  than 
the  corresponding  element  in  the  noise  matrix. 

Then, 

S  (uu)  =  S  (a>)  +  S  (oj)  =  S  (uu)  (59) 

—l  — s'  — n  — n 

and 

F  (uu)  =  s’l(uu)S  (uu)  S  ^(uu)  (60) 

—A  — n  — s  — n 

and 

G(uu)  =  Sn '(uu)  Sg(uo)  (61) 

The  performance  index  d2  becomes: 


{TrJ  (.)y«)Sa  W^W 

d3  *  2  — r - (62> 

{TrJ  S,-,(.)y.)Sj"(.)Ss(.)|} 
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Cancelling  the  common  term,  we  have: 


da 


1 

2 


Tr 


S  1  (uu)  S  (ou)  S  (id)S  (<d)  — 
J  n  —n  s  2tt 


(63) 


2.  Homogeneous  Noise, Threshold  Case 

In  this  case,  the  signal  spectrum  at  each  hydrophone  is  identical  and 
the  noise  spectra  are  also  identical .  Then  we  may  write 


S  (id)  =  S  (uu)  P(id) 

— s  s 

(64a) 

sn(uu)  =  sn(w)  g(uu) 

(64b) 

Here  P(d)  and  Q(uu)  represent  the  cross -spectral  terms. 
Then, 


and 


s,w 

SnW 


g'NaiJPOdQ'1^) 


(65) 


Ss^  -i 

G(id)  =  Q  (id)  P(id) 


(66) 


and 


d2  =  i 


s> 

S>> 


Tr  [q'l(«)P(«1)Q'l(«i)P(a»)]  ||  (67) 
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3 .  Single  Signal  Case 


We  now  turn  to  the  case  of  principal  interest  to  us.  Consider  the  array 
shown  in  Figure  7.  The  coordinate  of  each  hydrophone  is  denoted  by  a  vector,  r^. 
We  assume  that  the  signal  component  at  ^ach  hydrophone  is  identical  except  for  a 
time  delay. 


s(t  -  tx) 
s(t  -  t2) 


s(t)  = 


(68a) 


s(t  -  tn) 


ARRAY 


FIGURE  7  HYDROPHONE  ARRAY:  PLANE  WAVE  SIGNAL 
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A  simple  case  in  which  this  assumption  is  true  is  the  case  of  a  plane 
wave  whose  line  of  propagation  is  along  a  vector  a.  In  this  case 

a  •  r. 

T.  =  - -  (68b) 

1C 


We  will  consider  this  case  in  the  sequel  because  of  its  easy  physical  interpreta¬ 
tion.  The  modification  to  account  for  the  general  case  described  by  Eq.  (68a)  is 
obvious . 


Substituting  Eq.  (68b)  into  Eq.  (68a),  we  obtain 


s(t)  = 


(68c) 


Then, 


where 


S  (ui)  =  S  (uj)P(uu) 

S  S  "" 


(69a) 


+  )7S-t£t-r,] 


PKW  =  e 


(69b) 


In  this  case,  it  is  easier  to  perform  a  preliminary  operation  on  the 
inputs  so  that  the  outputs  of  the  hydrophones  are  in  time  synchronization.  This 
corresponds  to  "steering"  the  array  and  is  obtained  by  a  set  of  delay  lines,  as 
shown  in  Figure  8.  Observe  that  we  have  indicated  negative  delays;  these  are  ob¬ 
tained  physically  by  including  a  common  positive  delay.  The  noises  at  the  output 
are  a  function  of  the  stee^'.ng  direction. 
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and 


FIGURE  8  ARRAY  STEERING 


The  new  noise  spectral  matrix  follows  easily. 


ij 


S  (uu)  e 
— n. . 
ij 


,  .  U)  / 

+  j  -  (a.  •  r. 
c  -]  -j 


(70a) 


s' M  =  S  (uu) 


1  1  1 


(70b) 
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From  this  point  on  we  will  regard  the  steered  outputs  as  our  basic 
input.  For  notational  simplicity,  we  will  leave  off  the  primes. 


To  realize  the  filter  as  shown  in  Figure  3.  we  find  G04  as  given  by 
Eq.  (54).  Denote  the  elements  in  the  inverse  of  Si  04  as  SV  (4«  Then, 


GO  4  = 


— 

““ 

Si 11  (4  s il?(4  . 

„lN,  V 
•  Si  (id) 

s  (4 

9 

Ss(4  • 

•  Ss(,Ju) 

N1 

Si  (4 

s  (4 

S 

s  (4 

s 

S  (uu) 

s 

_ 

— 

_ 

(71a) 


or 


G(jx)  = 


N 

N 

N 

r 

L 

s?(4 

Z Silj(w)  •  • 

•I 

s\j(4 

j  =  i 

j  =  » 

j  =  i 

N 

I 

Sij  (4 

(71b) 

J  =  1 

N 

N 

r~- 

> 

L 

SiNj(4 

V 

l 

SfNj(4 

Lf-‘ 

J  =  1 

J 

The  output  is  an  N  x  1  matrix  >^04  which  is  defined  by  Eq.  (54b). 

Substituting  Eq.  (71b)  into  Eq.  (54b).  we  find  that  each  element  in  the 
matrix  is  identical.  The  elements  are  the  estimates  of 


}  (x)  =  S  04 

s  s 


N 

i 

L 


N 

£  SiiJ04  R.(j4 


i  =  i  j  =  i 


(72) 
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From  Figure  6,  we  see  that  we  want  to  form  the  dot  product  between 
z(uu)  and  £s(jju).  In  Figure  9(a)  we  show  the  receiver  for  this  particular  case 
using  matrix  notation.  In  Figure  9(b).  we  show  the  actual  operations.  Note  that 
each  element  in  ^>s(jai)  is  the  same.  Therefore,  we  can  add  the  elements  in  z(u) 
and  then  multiply  the  two  scalar  quantities.  The  resulting  receiver  is  shown  in 
Figure  10. 


(a) 


FIGURE  9  OPTIMUM  ARRAY  PROCESSOR:  SINGLE  SIGNAL  SOURCE 
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FIGURE  10  OPTIMUM  ARRAY  PROCESSOR:  SINGLE  SIGNAL 
SOURCE (ALTERNATE  FORM) 


The  difficulty  with  the  receiver  structure  in  Figure  10  is  that  there 
are  two  separate  combining  operations.  We  will  eliminate  this  difficulty  by 
proving  that  we  can  obtain  >s0i)  by  passing  z(j’*)  through  a  (scalar)  filter,  FGx). 
In  other  words,  we  want  to  prove: 

;sg*>  =  f(j*)  *0-*;)  c?3a) 

or 

N  N  N  N 

S  (*)  SiJ(ji)  R  0*>  =  F(jai)  /  y  Sij(ji;)  R  (j*)  P3b) 

S  i—  «_  J  L-.i-.Ti  } 

iri  J  =  i  i3*  Jsl 
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In  matrix  notation, 


SO.)  I  ^  O'-u)  RO'Ju)  =  _l^F(j'u)S  (j'Ju)  RG'x) 


(73c) 


where 


Equivalently,  we  must  prove: 


1  =  fl  1  .  .  .  1] 


(74) 


SO*)  1^  Si  0*)  ~  IT  F(j-ju)  S  0*) 
s  —  —  ~n 


(75) 


Post -multiplying  by  0  Ju).  we  have 


S(]J))1T  =  j.TF0r)s  'X0^)  Si(jju) 


(76) 


Using  Eq.  (70b)  and  the  definition  of  Si(tu),  we  obtain. 


SOx)  j_T  =  1TF0jj)S  Ox) 

s 


N 

1  slj 

L  n 

j  =  i 
N 

r  i!i 

/  s 
<->  n 

N 

7  sN) 

L.  n 
Jal 


N  N 

Vs’-j  ...  )  s1-* 

L  n  L  n 


i  =  i 


j  =  i 


(77) 


N 

7  S 


Nj 


♦  FOx)  lTl 
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This  Implies 


Ss(jj>) 


Defining 


N 

F(j*)S»  ; 

S  i _ . 


N 

;  s‘Jo-x)  +  F(jx) 

L*  n 


i=i  j=i 


N  N 

aV)  £  r  s‘JOx) 

4 n 
i=i j=i 

we  have: 


F(i)  = 


A  O  x)  S  O  x) 
s 


S  Ox)  +  A  Ox) 
s 


(78) 


(79) 


(80) 


The  resulting  realization  is  shown  tn  Figure  11.  Alternately,  we  see 
that  we  can  realize  the  receiver  in  the  filter -squarcr  form  as  shown  in  Figures 
1 2a  and  12b.  (The  purpose  of  the  modification  in  Figure  12b  will  be  apparent  in 
Chapter  IV.)  In  Figure  13.  we  show  the  resulting  structure  under  the  threshold 
assumption.  (Here,  we  assume  Ss(xi)  <  <  A  Ox)  for  frequencies  of  interest.) 

This  completes  our  discussion  of  the  general  receiver  structure  for 
the  detection  problem.  In  Chapters  V.  VI.  and  VII.  we  will  look  at  some  particu¬ 
lar  problems.  Before  doing  this,  we  will  demonstrate  how  a  similar  structure 
arises  when  a  waveform -estimation  approach  is  used. 


FKICRH  !  1  SIMPLIFIKD  FORM  OF  OPTIMUM  ARRAY  PROCFSSOR 
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(a) 


(b) 

FIGURE  12  FILTER  SQUARER  VERSION  OF  OPTIMUM  ARRAY  PROCESSOR 


l 


FIGURE  13  THRESHOLD  CASE 
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IV.  SIGNAL  WAVEFORM  ESTIMATION 


In  many  cases  detection  is  not  tne  only  question  of  interest.  For 
example,  we  may  want  to  process  r(t)  to  obtain  a  good  estimate  of  s(t)  to  use 
in  a  classification  problem. 

In  this  chapter,  we  discuss  two  possible  approaches  to  waveform 
estimation.  The  first  is  called  "distortionless''  filtering  because  the  processor 
is  designed  so  a not  to  distort  the  signal.  The  second  is  tne  classical  mini¬ 
mum  mean-square  error  approach.  After  deriving  them  separately,  we  show 
that  a  matrix  MMSE  filter  can  be  viewed  as  a  cascade  of  a  matrix  distortion¬ 
less  filter  and  a  sca!ar  MMSE  filter.  This  result,  which  is  due  to  Kelly, 
is  important  because  it  enables  us  to  perform  a  distortionless  combining  of  the 
hydrophone  outputs  and  observe  it  before  introducing  the  signal  distortion  caused 
by  an  MMSE  filter. 

Finally,  we  relate  the  waveform  estimation  problem  to  the  detection 
problem  and  show  that  tV  distortionless  combining  operation  :s  the  only  matrix 
operation  necessary  in  both  cases. 

In  this  chapter  we  will  assume  that  the  signal  is  a  plane  wave  whose 
direction  is  known.  As  discus.-..*.!  in  Section  B-3  of  Chapter  III,  we  assume  that 
.lie  delays  steer  the  array  on  the  target  have  already  been  inserted. 

rhus,  r  (t)  is  an  N-element  column  matrix  whose  elements  are:^ 

=  s(t)  +  n.(t)  (81) 


The  spectra)  matrix  is  S^uj),  where: 


CO 

S  (x)  =  j  E  [  n  (t)  n^(t  -  t)1  e  dT 

—  n  J  —  — 

-  » 


(82) 


with  this  mode!,  we  first  develop  the  idea  of  distortionless  filtering. 


tA  s  before,  the  primes  are  suppressed  for  notational  simplicity. 
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A.  DISTORTIONLESS  FILTERS* 

The  filter  of  interest  is  shown  in  Figure  14.  The  outputs  of  the  beam¬ 
forming  operation  are  given  by  Eq.  81 .  Each  output  contains  the  signal  s(t)  and 
a  noise  n j (t ) .  These  inputs  are  filtered  and  summed  to  give  a  single  output  y(t). 


O 

o 

o 


FIGURE  14  DISTORTIONLESS  COMBINING 


It  is  required  that,  in  the  absence  of  noise, 

y(t)  =  s(t) 


(83) 


for  any  signal  s(t). 

Under  this  constraint,  we  wish  to  minimize  the  variance  of  nc(t). 


tThe  idea  of  combining  multiple  inputs  in  a  statistically  oDtimum  manner  under 
the  constraint  of  no  signal  distortion  is  due  to  Darlington. (^)  An  interesting 
discussion  of  the  method  is  contained  in  Brown  and  Nilsson's  text. (20)  ft  was 
deri’.ed  independently  by  Levin^^Oas  a  minimum -variance  unbiased  estimator. 
A  simple  derivation  is  given  by  Kelly; 8)  our  discussion  follows  this  later  ref¬ 
erence  . 
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The  output  of  the  filter  is: 


Y(juu)  =  D  (ju))R(jtu)  =  £  D.O-Ju)  R  O-d) 


where  D(joi)  is  N  x  1  matrix. 


P(juo)  = 


D^uu) 

d2G'ju) 


dno^) 


The  constraint  of  no  distortion  implies: 


l  Di «»> " 


(86a) 


lT  D(jui)  =  DT0.)  1  =  1 


(86b) 


for  all  uu. 

The  variance  is: 


°’n  -  E["cft>]  *  J  S„  «">  f 
c  -  <*>  c 
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Now, 


S  (<jj)  =  E  D  0^)  n  Gw)  n  Gw)  DG<J>)  I 

n  L—  -c  ~c  —  _! 

c 


(88a) 


S  (uu)  =  D^Gw)  S  Gw)  D*Gw) 
n  —  — n  — 


(88b) 


We  want  to  minimize  c£c  subject  to  the  constraint  in  Eq.  (86b);  to  do 
this,  we  use  the  Lagrange  multiplier  technique: 


F  = 


|dtGw)  s^Gw)  d*Gw)  +  ^iTpOw)}^ 


(89a) 


Since  the  integrand  is  always  positive,  we  can  minimize  at  each 

frequency. 


_  A 
t  - 

UU 


pT0w)S  Gw)D*Gw)  +  XDTGujM 


(89b) 


obtain 


Differentiating  with  respect  to  DrGu),  the  rth  element  in  DGw),  we 


[0...1...]  SnG(i))D*Gw)  +  X  |0 _ 1 _ ]_i  =0 

r  =  1,2 _ N 

(90a) 


This  is  equivalent  to  a  single  matrix  equation. 


S  Gw)  D*Gw)  =  -  X  1 


(90b) 
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-1 

Pre -multiplying  by  Sn  Gw),  we  obtain 


D*Gw) 


-l 


XS  Gw)  1 


or 


r  -l  iT 

=  -  x  [sn  Gw) )  i 

or 

N 

D.Gw)  =  -  x  V  SljGw) 

j  L  n 

i  -  1 


(91a) 


(91b) 


(91c) 


The  value  of  X  is  obtained  from  the  constraint  equation  (86a): 

N  N  N 

\  D^w)  =  -  X  £  £  S*j  Gw)  =  1  (92) 

j  =  l  i=U  =  l 


so 


-  X 


N 


N 

c— '  < — i 

i-i  j*i 


-i  -l 


s1JGw) 

n 


=  A  Gw) 


(93) 


Therefore, 


N 

D.Gw)  =  AGo))  )  S^Gw) 

j  l  n 

i  =  1 


(94a) 
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The  distortionless  combiner  is  shown  in  Figure  15. 


FIGURE  15  OPTIMUM  COMBINER 


We  see  that  the  signal  z(jx)  is  identical  to  the  combiner  output  in 
Figure  11.  Similarly,  we  see  that  the  distortionless  waveform  Y(j'i))  also  ap¬ 
pears  in  Figures  12  and  13b. 

N  N 

l  l  s>>Rjw 

-  (94b) 

l  l  s> 

i=l j=l 


The  variance  using  distortionless  processing  follows  easily.  Substi¬ 
tuting  Eq.  (94)  into  Eq.  (91b)  and  the  result  into  Eq.  (88b),  we  have: 


S  (w)  =  /\(uu)lTS  (jw)S  (j^)  S  0^)1  Mw)  (95) 

n  -  — n  —  n  - n  — 

c 

But 


1T  S  lGw)'l  =  (A(jw)) 

-  — n  — 


(96) 
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Therefore, 


S  (ju)  =  A  Gw)  (97) 

n 

c 


and 


j  S  *  °c  <98> 

-  cc 

Several  observations  now  follow  easily: 

(1)  Looking  at  Figure  11,  we  see  that  s(t)  is  obtained  by 
operating  on  the  distortionless  output  y(t)  by  a  filter 


S 

s 


S  (<d)  +  A  (uu) 
s 


(99a) 


This  is  precisely  the  optimum  unrealizable  Wiener 
filter  for  estimating  s(t)  using  a  minimum  mean- 
square  error  criterion  when  the  input  is  y(t).  (Since 
the  "noise"  spectrum  is  A (d)) .  The  mean-square 
error  is 


mmse 


®  AH  S  (id) 

S 

J  s  (ui)  +  A  (id) 

®  s 


dip 

2tt 


(99b) 


(2)  All  of  the  processing  up  to  y(t)  in  Figures  12,  13b, 
and  16  is  independent  of  the  signal.  The  effect  is  a 
noise -reduction  due  to  combining.  In  Figure  16  we 
re-draw  the  receivers  in  Figure  12  and  13b  to  empha¬ 
size  the  distortionless  signal. 

Before  leaving  waveform  estimation,  we  will  consider  the  problem 
from  a  different  viewpoint. 
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B,  MAXIMUM  LIKELIHOOD  ESTIMATES  AND  GENERALIZED  LIKELIHOOD 
RATIO  TESTS 


The  derivation  of  the  optimum  distortionless  filter  started  with  a 
rather  arbitrary  linear  filter  restriction.  One  can  also  approach  the  problem 
from  the  standpoint  of  a  maximum  likelihood  estimate  of  s(t)J 

The  likelihood  function  of  r(t).  given  s(t),  is 


In  A  (r(t):s(t))  =  +  | 


*%)  IT  S' \u)  R(i))~ 

—  n  —  ZTT 


1 

+  2 


R*T(uu)S  l(u)  l 
—  n  ~  2tt 


(100a) 


CD 

*  J  1  (uu)  i#(uu) 

-  00 


where  we  have  assumed  an  infinite  observation  interval  and  stationary  processes. 

A 

The  maximum  likelihood  estimate  of  $(ju)  is  simply  the  function  j*m^(ju) 
that  maximizes  the  likelihood  function.  To  find  it,  we  set 

Hu)  =  J  ,(x)  +  e  (uu)  ,  (100b) 

mt  c 

differentiate  In  A  with  respect  to  e,  and  require  the  result  to  equal  zero  for  e  =  0 
and  all  *  (i ).  This  gives 


*  djj  r  T  -1  T  - 1 

(x-)5§  {lS  (x)  Kh)  -  I  S  (x)l  /  .(x’)f 
.  •:  2  t  —  - n  —  -  n  -  mt  j 

(100c) 

-  co 

tThe  r-  ,-uits  through  Eq.  lOOe  were  previously  derived  by  Kelly  and  Levin. ^2) 
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The  two  bracketed  terms  are  conjugates.  Since  «*e(i})  is  arbitrary,  the  term 
inside  the  bracket  must  equal  zero  for  all  ju.  Therefore, 


«i\W 


ITS  (uu)  R(ju) 


(lOOd) 


Writing  Eq.  (100d)out,  we  have: 


i 


N  N 

L  L  s>>Vr) 

LliJji _ 

N  N 


L  i 


s> 


i  =  i  j  =  i 


(100e) 


A 

Looking  at  Eq.  (94b).  we  see  that  >  ,(jju)  is  identical  to  Y(jx). 

mt 

Observe  that  there  was  no  a  priori  restriction  to  a  linear  processor. 

The  detection  analog  to  the  maximum  likelihood  estimate  is  the  gen¬ 
eralized  likelihood  ratio  test  (e.g. ,  Davenport  and  Root,  Chap.  14^)  or  Van  Trees, 
Section  2.5.^)) 

The  generalized  likelihood  ratio  test  is 

max  p  jjr(t) | 

ft  (r(t»  =  - - -  >  v  (101a) 

8  p  LI(t)l  H0j 


The  numerator  is  simply  p  [r(t)l  Hjl  evaluated  at  s(t)  =  snij(t) 
Substituting  Eq.  (lOOe)  into  Eq.  (lOOa)gives  {n  i\g(r(t)). 


49 


3rthur  3I.ittlr.Jnr 


'n  .'  g(r(t)) 


At)  1T  S*l(t)  R(4^ 
m  ’  -  -n  -  2” 


(101  h) 


t 

+  2 


(vVe  incorporated  the  last  term  into  the  threshold.) 
Now. 


>  ,(x)  =  A(x)  lTS*‘(-j»)  R(i-) 
m  >  —  -n  - 


=  Ml)  RT(jl)S  (4  1 
—  — n  — 


(101c) 


=  A(4  R  (ju)S  ‘(4  1 


Substituting  Eq.  (101c)  into  Eq.  (101b).  we  obtain 


1  "  -•  djL 
In  A g  (r (t ))  +  j  A(x)  2" 


1  "  -  .  .  dx 

2  .  A(x)  *m '(l) 2  2~ 


( 1 0  kl ) 


(Recall  A(i)  is  real) 


SO 
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Therefore. 


<  n  Ag(_r(t))  = 

-  ® 


(x) 


tli 


(I  Ole) 


Using  Parseval's  theorem,  we  obtain  the  receiver  for  the  generalized 
likelihood  ratio  test  that  is  shown  in  Figure  17. 


MATRIX 

A 

T/J 

FILTH 

1 

SOUARER 

r  * 

j 

-Vi 

LINEAR  FILTER 

FIGURE  17  RECEIVER  TO  IMPLEMENT  GENERALIZED 
LIKELIHOOD  RATIO  TEST 


C.  SUMMARY 


In  this  chapter,  we  have  related  the  ideas  of  distortionless  filtering 
and  maximum  likelihood  estimation.  Criteria  of  this  type  are  appropriate  when 
the  signal  is  a  non-random,  but  unknown,  waveform.  Further,  we  found  that  the 
generalized  Iikelih<xxl  ratio  test  is  closely  related  to  this  estimation  procedure. 

Combining  these  result-  with  those  of  Chapter  HI  we  obtain  the  com 
posite  receiver  shown  in  Figure  IS. 


Si 
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BAYES  OR 

NEYMAN-  PEARSON 
TEST _ _ 


FIGURE  18  CRITERION  INVARIANT  RECEIVER 
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V.  HOMOGENEOUS  NOISE:  ARRAY  GAIN 


A  special  case  of  interest  is  when  the  noise  is  homogeneous. 
Specifically, 

S  (jw)  =  S  (ju»)  =  S  (j'4  (102a) 

n,,  n..  n 

11  n 

and 

Sn  ^  =  Sn<J,J,>  PijOx)  (102b) 

n.j  n  ij 


Let 


Then, 


£(juu)  be  a  normalized  cross -spectral  matrix  with 
elements  Pij(jjj) 

S^Guu)  =  (SnGui))  1  P  (jiu) 


(102c) 


One  realization  of  the  resulting  receiver  is  shown  in  Figure  19. 

In  this  case,  it  is  easy  to  evaluate  the  effect  of  the  array  on  the  per¬ 
formance  of  the  system  with  respect  to  the  following  three  criteria: 

(1)  Detection  performance  index  (Eq.  (67)) 

(2)  Distortionless  signal:  minimum  noise  variance 
(Eq.  (98)) 

(3)  Minimum -mean-square  filtering  error  (Eq.  (99b)) 

We  shall  consider  these  three  cases  in  order .  At  a  single  hydrophone: 


uu 

'TI J 


Ssw 

SnM 


duo 

2tt 


(103) 


This  is  simply  the  scalar  version  of  Eq.  (67) 
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To  evaluate  d2,  we  look  at  the  signal  spectrum  of  y(t)  in  Figure  I9.t 
The  signal  spectrum  is  Ss(uu),  and  the  noise  spectrum  is  A(w).  Therefore, 


but 


d2 


duo 

2rr 


A  1(uu) 


N  N 

=  1  l 

i=i  j=i 


p1J(jw) 


s  (id) 
n 


A  A(uj) 

-  snw 


Therefore, 


(104) 


(105) 


d2 


s^uo 

Sn2(uj) 


A2(uj) 


dijj 

2tt 


(106) 


We  see  that  the  effect  of  the  array  is  contained  completely  in  the  func¬ 
tion  A(uu).  This  function  is  commonly  referred  to  as  the  array  gain 


s  n  N  N 

A(uu)  =  -773  =  Y  £  Pij  («) 


A(uj) 


(107) 


i=i j=i 


For  independent  noises,  the  array  gain  is  simply  N. 

We  recall  that  the  elements  p^(ui)  incorporate  the  effects  of  array 
steering.  In  terms  of  the  unsteered  noise  matrix, 

N  N  .  w  r  -I 

-j  —  ol  •  l.r.  -  r.  J  , . 

c  -  -1  -j  J),"A  (108) 


Ato)  =  11  e 


*:"» 


1  J 


tWe  could  also  work  directly  with  Eq.  (67)  (see  Appendix  C). 
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Thus,  unless  the  noise  matrix  is  diagonal,  the  array  gain  is  a  function  ui  the 
steering  angle . 


The  array  gain  also  arises  naturally  in  the  other  cases . 

For  distortionless  filteiing,  the  variance  at  a  single  hydrophone  is: 


-  J  S„<”>  R 


(109) 


At  the  combiner  output. 


■  I 


- 


S  (id)  , 
n  dii) 

A(id)  2tt 


(110) 


For  minimum -mean- square  error  filtering,  the  error  at  a  single 
hydrophone  is: 


J 

-  t 


S  (id)  S  (uu) 
s  n 

S>)  +  S  ((«) 
s  n 


dtp 

2tt 


At  the  final  output, 


S  (u ))  A  (id) 
s _ 

S  (id)  +  A  (id) 
s 


dip 

2tt 


dU' 

2tt 


(HD 


(112) 


Therefore,  for  any  of  the  three  purposes,  the  array  gain  completely 
characterizes  the  effect  of  the  array.  It  is  dependent  on  frequency  and  steering 
angle,  and  its  effect  is  to  reduce  the  noise. 

In  the  next  section  we  derive  various  noise  models  and  find  the  array 
gain  for  some  interesting  configurations. 
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VI.  DISTRIBUTED  NOISE  FIELDS 


In  Chapters  VI  and  VII  we  apply  the  preceding  results  to  simple 
examples.  In  this  chapter,  we  shall  consider  noise  fields  that  are  distributed 
over  a  region  in  space .  * 


A.  HYDROPHONE  NOISE  ONLY 


Example  1 


In  this  case,  the  only  noise  is  hydrophone  noise,  which  is  assumed 
to  be  independent  from  one  hydrophone  to  the  next.  For  simplicity,  we  first 
assume  that  the  noise  spectra  are  identical .  Thus , 


S  . .  (id)  =  S  (id)  6..  (113) 

n.ij  n  ij 

In  this  particular  case,  it  is  clear  that  the  spectra  after  steering 
are  identical . 

Thus, 


(114) 


Clearly, 

i>  =  sTT>  i  <U5> 

n  ' 

Combining  consists  of  a  simple  summing  operation. 

The  array  gain  equals  the  number  of  elements  in  the  array: 


A  (uj)  =  N  (116) 

Since  all  of  the  noise  is  assumed  to  be  hydrophone  noise,  the  spacing  and 
location  of  the  elements  are  unimportant. 

The  modification  to  include  different  noise  levels  is  straightforward. 

*As  pointed  out  in  the  references.the  results  in  this  chapter  are  due  to  Dr.  E.J. 
Kelly,  Jr.  His  contribution  is  important  to  the  over -all  unity  of  the  report,  and 
we  are  happy  to  acknowledge  it . 
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B.  ISOTROPIC  NOISE  ONLY 


Several  models  of  noise  geometry  lead  to  what  is  called  an  isotropic 
noise  field.  Two  of  these  are: 

(1)  The  noise  is  assumed  to  consist  of  uncorrelated 
plane  waves  with  identical  statistics  coming  from 
all  directions  (e.g.,  Marsh(23)). 

(2)  The  noise  is  generated  by  uncorrelated  noise 
sources  uniformly  distributed  on  the  surface  of 
a  large  sphere  (e.g.,  Faran  and  Hills (24) 

If  the  noise  is  assumed  to  consist  of  a  single  frequency,  then  it  is 
easy  to  3how  that 

P .  .  ( T )  =  Sln  ldU.  COS  2TTfQT  (117) 

1J  kd. . 

iJ 

where 

2nfo 

k  =  - -  is  the  wave  number 

c 

c  is  the  velocity  of  sound 

f  is  the  frequency 

and 

d. .  is  the  distance  between  elements 

D 

Due  to  the  spherical  distribution  assumption,  the  correlation  func¬ 
tion  (before  steering)  does  not  depend  on  the  orientation  between  the  two  ele¬ 
ments  . 


A  simple  extension  of  the  single -frequency  case  is  when  the  noise 
spectrum  at  each  point  on  the  sphere  is  the  same  (say  Sn  (w)  ).  Then,  it 


follows  easily  that 


0. .  (jn) 

lj 


d.. 

in  uu  — ^ 
c 

c 


S  (U!) 
n 


(118) 


The  correlation  function  may  be  found,  but  it  is  not  necessary. 
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Kxample  2 


To  illustrate  the  isotropic  noise  case,  we  consider  a  simple  ex 
ample  due  to  Bryn(^),  The  array  is  shown  in  Figure  20. 


SIGNAL 

« . - 


I 

1 

FIGURE  20  UNIFORMLY  SPACED  FIVE-ELEMENT  ARRAY 

It  consists  of  tive  elements  arranged  in  a  line.  The  signal  direction  is  per¬ 
pendicular  to  the  array.  In  this  case,  no  delays  ar®  necessary  for  steering. 
The  normalized  spectral  matrix  is: 


where 


p(jw)  - 


,  sin  u 
u 

sin  2u 

2u 

sin  3u 

3u 

sin  4u 
4u 

i 

sin  u 

u 

sin  2u 

2u 

sin  3u 
3u 

1 

sin  u 

u 

sin  2u 

2u 

(symmetrical) 

1 

sin  u 

1 

2nfd  _  d 

-  =  2tt  — 

c  X 


(119a) 


(119b) 
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The  array  gain,  A0  (f),  as  a  function  of  u  is  shown  in  ligure  21 . 
The  gain  obtained  by  a  straight  sum  is  also  shown  in  the  figure.  We  see  that 
for  spacings  greater  than  half  a  wavelength  the  two  are  esseniOlly  identical . 


(from  Iryn;  rt f.  13) 


HGURfc  21  ARRAY  GAINS 


In  this  model  the  gains  correspond  to  those  obtained  by  Pritchard'25) 
using  a  signal -to -noise  ratio  criterion.  It  turns  out  that  the  gains  are  'arge 
and  opposite,  inis  type  oi  array  is  commonly  referred  to  as  a  "super-gain" 
array:  as  one  would  expect,  such  arrays  are  sensitive  to  variations  in  the 
gains  and  therefore  are  difficult  to  realize  in  practice,  if  one  considers  more 
general  arrays  (the  elements  not  necessarily  ir  a  line),  then  one  can  get  a 
singular  detection  problem  for  the  isotropic  noise  model  (see  Gaarder'^b)  ). 
One  can  eliminate  this  sensitivity  and  the  possibility  of  singularity'  by  keeping 
the  hydrophone  noise  it  a  non -zero  ievel . 

Other  array  gain  calculations  have  been  made  (e.g.,  B.  Cron  and 
C.  Becker<  7)  ). 

We  now  turn  to  a  second  category  of  noise. 
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VII.  DIKi  C  1 1  >N  A  L  M  US  I  M)I  K<  IS 


I- requcntlv ,  the  noise  th.it  we  want  to  combat  has  a  strong  directional 
character.  The  limiting  form  of  this  category  is  a  noise  plane  wave  corning 
from  a  single  direction.  We  investigate  this  Lmiting  case  under  the  assump¬ 
tion  that  ihe  direction  is  known  exactly.  If  the  direction  were  not  known,  the 
receiver  would  have  to  measure  it.  The  results  for  the  known  direction  pro¬ 
vide  a  bound  on  how  well  a  receiver  incorporating  measurement  could  do. 

A  simple  model  of  directional  interference  is  shown  in  Figure  22. 


l 

i 

♦ 

I 

i 

l 

: 


FIGURE  22  DIRECTIONAL  NOISE  MODEL 


The  noise  is  a  stationary  random  process  from  a  known  angle  s. 
For  simplicity,  we  assume  that  the  signal  direction  is  broadside.  The  arra\ 
elements  are  uniformly  spaced  in  a  line.  If  we  label  the  noise  voltage  at  the 
center  element  n(t) ,  then  lor  an  array  with  2N  1  elements. 
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where 


n  (t  +NT) 

• 

n  (t  +  2  t  ) 
n(t+T) 
n  (t) 

n  (t  -  t) 

• 

n  (t  -  Nt  ) 

t  =  d/c  sin0 


(120) 


The  noise  cross -spectral  matrix  is  easy  to  find.  For  a  three - 
element  array, 


r. 


e 


j  OUT 


e 


-juu2T 


(01)  =  Sd(uu) 


+  jUJT 

e 

+jll)2T 


1 


(121) 


e+juj-r 


1 


It  is  clear  that  this  matrix  is  singular.  For  any  value  of  9  except 
zero  we  can  achieve  perfect  detection.  Under  ,  for  example, 

rl  (t_T)  '  r2  (t)  =  s(t_T)  -  s(l)  (122a) 

and'inderHg,  ^  (t-T)  -  r2  (t)  =  0  (122b) 

which  gives  a  perfect  detection  capability . 

Therefore,  to  make  the  pure  directional  noise  problem  meaningful 
we  must  include  a  white  noise  component.  Then, 

N0 

sn  to)  =  I  +  s  (U!)  (123) 

— n  2  —  n , 
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For  the  three -element  case, 


S  (ou) 


=  (t+V-">) 


where 


x  = 


1 

-jX’T 

xe 

xe 

-ju)2r" 

-Hui’T 

xe 

1 

xe 

-jUJT 

+j  2  xT 
xe  3 

+jjUT 

xe 

1 

- 

Sd<*> 

(124) 


N0 

T  +  V*> 


(125) 


The  array  gain  is  the  sum  of  all  the  elen^  ats  in  the  inverse  of  the 
matrix  in  Eq.  (124).  Taking  the  inverse  and  summing  the  elements,  we  obtain: 


A  (jut) 


3  (1  -  x“)  +  2x  (x  -1)  [2  cos  out  +  cos  2  x't] 

!  +  2x3  -  3x2 


(126) 


As  we  would  expect,  ior  x  =  0 


A  (j  j.)  =  3 

To  indicate  the  behavior,  we  have  plotted  A  (j  u)  for  two  values 
of  uut  in  Figure  23, 

For  any  given  element  spacing,  these  can  be  translated  to  a  particular 
value  of  9  . 

For  example,  if  d  -  \  /2  ,  then 

'jut  =  n  sin  9  (127) 

and  the  two  curves  shown  represent  9  =  0  and  30°,  respectively.  We  have  also 
indicated  the  array  gain,  Ar(jju)  for  a  conventional  array  which  sums  the  outputs. 

The  procedure  for  finding  the  array  gain  for  larger  arrays  or  multiple 
directional  noise  sources  is  conceptually  straightforward,  and  the  actual  com¬ 
putation  can  be  done  numerically. 
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FIGURE  23  ARRAY  GAIN  VS 
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VIII.  CONCLUSIONS 


*■  / 

The  primary  purpose  of  this  report  has-been  to  fit  various  tech¬ 
niques  and  criteria  related  to  optimum  array  processing  into  a  unified 
theory.  Since  the  models  one-twes  are  only  approximations  to  the  actual 
physical  situation,  it  is  important  to  understand  how  various  assumptions 
affect  the  optimum  receiver  structure . 


For  an  interesting  class  of  criteria  and  signal  models  we  found 
that  the  optimum  receiver  consisted  of  a  set  of  delays  to  steer  the  array 
followed  by  a  combining  operation  which  depended  only  on  the  noise  covari 
ance  matrix.  The  output  of  this  combiner  is  a  single  waveform,  which  is 
then  processed  depending  on  the  criterion  and  signal  model . 
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APPENDIX  A 


DERIVATION  OF  EQUATION  21 


Using  Eq.  (14)  and  Eq.  (20),  we  can  write: 

hA(t,U)  =  -^(t.u)  +  Q0(t,u)  (A.  1) 

Pre -multiply  by  K^z.t),  integrate  with  respect  to  t,  and  use  Eq.  (9): 

Tt 

J  K1fe,t)h(t,u)dt  -  -1.6<z  -u)  +  J  K1(z,t)g0(t,u)dt  (A. 2) 
Ti  Ti  ’ 


Post-multiply  by  Kq(u,x),  integrate  with  respect  to  u,  and  use  Eq.  (9): 


j  lS1<z't)]lA<t*u)K0(u  z)dt  du  =  -  KQ(z,x)  +  K^z.x)  (A. 3) 

T. 

1 


Observing  that 


K^z.x)  -  Kq(z,x)  =  Kg(z,x) 


(A. 4) 


we  obtain  Eq.  (21) . 
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APPENDIX  B 


EVALUATION  OF  THRESHOLD 


We  want  to  evaluate  the  sum, 


(B.l) 


Using  Eq.  (16),  we  see  that 


(B.2) 


where  the  notation  2/Nq  emphasizes  the  dependence  of  h(  •  )  on  the  noise  level. 


One  can  show  that  error  matrix 


§  (t,t  :  z)  =  7=  i  [  h  (t,t  :  z) 

-C-  T  Z  J  -c- 


dt 


so 


(B.4) 


(B.5) 
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APPENDIX  C 


EVALUATION  OF  d2  BY  USE  OF  EQUATION  67 


In  Eq.  (104),  we  derived  d2  using  the  combined  signal.  Here,  we  use 
Eq.  (67)  directly. 


«  S» 


d2  =  |  J  Tr^'o^PO^^'O^PO^]^  (C.D 

-  co  n ' 


We  want  to  show  that 


Tr  |*Q  '(jw)  P(juu)  Q ^(juu)  P(jw)j  =  A2(u>) 

n  r  n  n 

<-  V  U. 

■  L  L  p  0uj) 


N  N 

I  I 

[i=l  j-l 


k=l  1=1 


(C.2) 


This  is  a  straightforward  exercise  in  matrix  manipulation. 
Recall  that,  after  steering, 


P(w)  = 


1  1  1  .  .  .  1 

1 

i  1  1  .  .  .  1 


and 


Q  1(jw)  = 


pll(w)  p13(w)  p13(u>)  PlN(u>) 


Ni.  , 
P  (t>) 


NN.  . 
P  (o) 


(C.3) 


(C  .4) 
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Therefore 


Q  *(jw)  P(jw)  = 


N  N 

^  piJ(u»)  £  pxW 
j=l  j=l 

N 

P2%) 

j  =  1 

N 

^  PNJ(.) 

j  =  l 


•  •  • 


N 

l  PNj(^) 
j  =  l 


(C5) 


Squaring  this,  we  obtain 


|q  Gw)  P0<jj)J  “ 


N  N  N 

I  Iplj  I’1 

k=  I  j  =  l  1  =  1 


kl 


X 


N 


I  I  p3j  I  P“ 

k= 1 j  =  1  l=i 


N  N  N 

I  I  p3j  I 

k=  1  j  =  l  l-  1 


k{. 


__J - 


(C.6) 


The  trace  is  simply 


N  N  N  N 


I  1  ;>lJ  Zp“-A"« 

i-l  k«i  j-1  t-l 


(C.7) 


which  is  the  deaired  result. 
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We  also  observe  that 


N  N 

I  »1J 

i=l  j=l 

so  that,  for  this  case, 


Tr  [Q'10  «)P0'JJ)]  =  £ 


(C  .8) 


{Tr  P(juo)3}  =  Tr  PQuu)}  j 


(C.9) 


The  expression  of  the  right-hand  side  of  Eq.  (C.9)  is  identical  to 
Bryn's  result. 
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GLOSSARY 


A'uu)  optimum  array  gain 

Ac(oi)  conventional  array  gain 

a  unit  vector  denoting  propagation 

direction 


c 


velocity  of  propagation 


d. .  inter -element  distance 

U 

d2  output  signal -to -noise  ratio 

p(ijy)  transfer  matrix  of  distortionless 

combiner 

6(t)  delta  function 


E  [•]  expectation  operation 

E.  [•  ]  expectation  assuming  i1*5  hypothesis 

is  true 

r\  threshold 


F(F  furctions  used  in  minimization 

f(t)  known  signal 

g(ui)  transfer  function  of  matrix  filter 

Y  threshold 


H0’  H1 


matrix  filters 

coefficient  in  series  expansion 
hypotheses 
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l_ 

identity  matrix 

1 

A 

unity  matrix 

k 

wave  number 

So-  Si-  5s’  S» 

covariance  matrices 

k(z,t) 

matrix  filter 

L 

sufficient  statistic 

A(r<t» 

likelihood  ratio 

v> 

generalized  likelihood  ratio 

X 

Lagrange  multiplier 

M«»> 

filter  function 

s 

scalar  eigenvalues 

m^t) 

mean -value  vector  on  Hj 

n(t) 

vector  noise  process 

nc(.) 

colored  noise  component 

nc(0 

combined  noise  component 

OJ 

frequency 

P(vJU) 

matrix  denoting  cross -spectra 

^(t) 

vector  eigenfunction 

inverse  matrix  kernel 

Q(w) 

matrix  denoting  cross -spectra 
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P(w) 

r(t) 

r.(t) 


normalized  cross -spectral  matrix 

received  vector  waveform 

components  of  received  vector 
waveform 


r. 

-l 


three-dimensional  position  vector 


R(w) 


Fourier  transform  of  r(t) 


s(t) 

si<(> 

SR(t) 

S(D 
s(t:  ii) 
S(a) 

Mw) 

o2 

c 

o2 

"mmse 


vector  signal 
illuminating  signal 

reflected  signal 

estimate  of  vector  signal 
random  signal 
spectral  matrix 
Fourier  transform  of  s(t) 
variance  of  combined  noise 

variance 

minimum -mean-square  error 


T 

1 


T  = 


e 

T. 

l 


T 

T 


f 


u 


delay 

argument  of  covariance  function 
phase  of  reflected  signal 

angle  of  directional  noise 
initial  observation  time 

final  observation  time 

length  of  observation  interval 


2  ^tL 

(  =  -T-— )  normalized  .ariable 
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VR 

amplitude  of  reflected  signal 

w(t) 

white  noise  waveform 

y(t) 

distortionless  output 

z(t) 

scalar  output  of  combiner 

matrix  notation: 

Tr  [  ] 

trace 

mt 

transpose 

l-  r1 

inverse 

* 

L*  J 

conjugate 

sij 

ij  element  in  inverse 
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